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Enlarged symmetry characterized by the group SU (4) can be realized in isolated semiconducting 
quantum dots. A Hubbard model then describes a pillar array of coupled dots and at half-filling 
the system can be mapped onto an S'?7(4) spin chain. The physics of these new structures is 
rich as novel phases are attainable. The spins spontaneously dimerize and this state is robust to 
' perturbations which break 5'(7(4) symmetry. We propose ways to experimentally verify the existence 

, of the dimerized phase. 
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^ ■ I. INTRODUCTION 

Qj ' Quantum dot arrays are a new arena for the study of strongly correlated electrons and the persistence of quantum 
coherence. Physical properties of a single semiconducting dot as well as tunneling between dots can be controlled over 
. ■ a wide range - a luxury not available to us in ordinary condensed materials. Recent advances in nanofabrication tech- 
I niques offer the possibility of constructing artificial structures so small that the electronic level spacing is comparable 
to the charging energy. As a consequence, these structures can exhibit enlarged continuous symmetries not normally 
I found in nature. In this paper we determine conditions under which a pillar of coupled semiconducting quantum dots 
: realizes the group 5'C/(4) as a good symmetry and show that the SU{4:) spins spontaneously dimerize - a phase of 
Ch ■ matter that would be difficult to attain with the smaller SU{2) symmetry of electrons in generic quantum dots. 

Continuous symmetries are ubiquitous in physics. Rotational invariance characterized by the group 0(3) permits 
the classification of atomic orbitals via integer angular momentum quantum numbersH. Spinning particles, such as 
electrons, are described by representations of the group SU{2). Approximate. 5[/(3) isospin symmetry of hadrons has 
i its origin in the light masses of the up, down and strange flavors of quarksa. Unlike the case of electrons for which 
SU{2) symmetry is exact, quarks can be described by SU(3) only approximately since the masses of the quarks are 
not exactly equal. Nevertheless, the approximate 5[/(3) symmetry is useful for classifying hadrons. We show how an 
' approximate S'i7(4) symmetry can be realized in quantum dot structures, and we exploit its properties to describe 
, the novel dimerized phase that should emerge in these structures under certain conditions. 

' Consider a potential well with N degenerate eigenstates. Taking electron spin into account there are a total of 
' 2 X N degenerate states, and if all of these states are equivalent, in a sense made precise below, we can think of them 
' as realizing the fundamental representation of the SU{2N) group. In other words, electrons placed in the shell can 
be considered as having 2N different, but equivalent, flavors instead of just the ordinary two flavors of spin up and 
down. It is important to note that SU{2N) symmetry is not equivalent, in general, to the higher-spin representations 
of the usual SU{2) group familiar from the quantum theory of angular momentum. Rather, for iV > 1, SU{2N) is a 
different and larger symmetry. 
' Ordinary atomic orbitals might seem like a good candidate but, for real atoms, the enlarged symmetry is broken 
Q , down to the usual symmetry by electron-electron interactions which lift the degeneracy. However, as Stafford 

O ' and Das Sarma noticedtl, semiconducting quantum dots offer the possibility of realizing enlarged symmetries. Quantum 
^ [ dots can be thought of as artificial atomsQEl with tunable parameters. To be precise, the electron mass is replaced 
' by the smaller band mass uie — > mi,, and the Bohr radius as = / [mf.e^) is replaced hy a*^ = e (me/mt,) as- In 
/\t ' GaAs, rrib w 0.067me, the dielectric constant e ~ 13, and « lOOA which is two orders of magnitude laifger than its 
fundamental value. Electrons in a quantum dot are confined in a non-singular potential often describedu as a short 
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square well in the z-direction and a simple parabola in the x ~ y plane, V{x,y) = ^ rUb uja^{x'^ -\- y^), though our 
results do not depend on the detailed form of the potential as long as it has cylindrical symmetry. For the lowest 
mode in the z-direction, the resulting harmonic oscillator eigenenergies are -En,;^ — 7iwo(2n -I- jZ^I -I- 1) where n and Iz 
are respectively the radial and angular momentum quantum numbers. 
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II. CONSTRUCTING AND COUPLING SU{2N) QUANTUM DOTS 



We-|propose a one-dimensional array of rotationally symmetric III-V semiconducting quantum dots arranged in a 
pillarQ and show how approximate SU{4:) syisnmetry can be realized in the structure. Given sufficient control over 
dot diameters and gate positions and biaseaS the lowest s-level (n = Iz = 0) of each dot can be completely filled 
(with two electrons) , rand the next higher, four- fold degenerate p-level with n = and Iz = ±1, can be half- filled with 
two valence electronsO, as shown in Fig. |^. Apart from configuration splitting (discussed below), the two p-electrons 
realize a self-conjugate (particle-hole symmetric) representation of SU{'i). The dimension of this representation is 
6 which corresponds to the six distinct ways the two electrons can be placed in the four available p-stateaJ. Once 
electron tunneling between the dots is turned on, there are four energy scales in the problem: the gross level spacing 
AE — Kluq in each well, the on-site Coulomb repulsion energy U which represents the energy cost to add an additional 
electron to the dot, the energy splitting between the six different p-level configurations AU, and the tight-binding 
electron hopping amplitude t > between states in adjacent wells, see Fig. |l| (B). 

(A) 

(B) 




FIG. 1. (A) Individual quantum dot and proposed pillar array of quantum dots with approximate 5(7(4) symmetry. By 
adjusting the bias, the lowest s-level in each dot is filled completely and the first p-level is half-filled with two electrons. (B) 
Energy diagram of two adjacent dots from the array. AU is the energy difference between the highest state of total orbital 
angular momentum \Lz \ = 2 (shown in the right) dot and the lowest state with = (depicted in the left dot). An electron 
in a p-level can temporary hop into an empty level on an adjacent dot and then hop back. This virtual exchange process lowers 
its energy by of order J = At^ /U. If !7 S> AU and J 3> AU, all 6 configurations on each dot participate equally in the exchange 
and the array realizes approximate SU (4) symmetry. 

The advantage of proposed pillar array Fig. |^ (A) is clear: Conservation of the electron's orbital angular momentum 
around the z-axis, a consequence of the cylindrical symmetry of the confining potential, guarantees that transitions 
between different angular momentum states in adjacent dots, which would break the flavor symmetry, are forbidden. 
The crucial conditions are that the energy gain J due to electron exchange between the dots greatly exceed AU, and 
also that the four flavors of electrons participate in the exchange on an equal footing. To second order in perturbation 
theoryli3 J — At^ /U and thus we require: 

— > A[/ and L/ > AC/ . (1) 

Another inequality ensures that only the p-electrons play an active role in the low energy physics: 

AE->t . (2) 

To estimate the size of AU for the p-states in a quantum dot we first note that spin-orbit coupling is negligibleS. Thus, 
AU is due almost entirely to the dependence of the electron-electron Coulomb interaction on the shell configuration 
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as described by Hund's rulest3, which have been shown both experimentallyQ and theoreticallyc^lEJ to be directly 
appUcable to semiconducting dots. The six configurations break up as: 6 ^ 3 © 1 © 2. The triply degenerate state of 
total orbital angular momentum = and total spin S" = 1 is lowest in energy, the intermediate non-degenerate state 
has Lz — and S = 0, and the two-fold degenerate highest level has \Lz\ = 2 and S = 0. As both AU and U scale 
as d^^ with dot size d, we may introduce T = /S.U/U, where T depends only the shape of the dot and the confining 
potential. Lowest order direct and exchange interaction integrals allow us to estimate that T ranges from 0.5 for thin, 
quasi-two-dimensional dots to F = 0.2 for thick dots. Because the Coulomb interaction is long-ranged, these numbers 
are nearly independent of the confining potential; indeed, T k, 0.2 also holds for real (nearly spherical) atoms. For an 
electron in a potential well of characteristic size d we have: AE h? /(rubd'^), U « e'^/{ed), and thus AE/U w a*g/d. 
Symmetry breaking effects due to the electron-electron interaction are therefore minimal in sufficiently small dots. To 
satisfy Eqs. and (|^) with F = 0.2, simple algebra shows that we require {AE/U)'r-S> 1 which is in fact satisfied by 
small dots. For example, InAs/GaAs {a*g « 300A) quantum dots have been madella which have d « 200A, t/ w 18 
meV, AE w 50 meV, and {AE/UY ^-p- Adjusting the array spacing and the thickness of the insulating barriers the 
hopping amplitude may be increasedlla io t = 0.2 AE. It then follows that J k, meV. For dots which arc not too 
thin, F — 0.2, AU ~ 4 meV, and the crucial inequalities Eqs. (|l|) are satisfied. In contrast to these artificial atoms, 
there are no real atoms for which both ineaualities Eq. and Eqs. (||) hold because AE sa U and t -^U . K typical 
example is a copper-oxide antiferromagnetO with J « 0.13 eV, U « 10.5 eV and hence AU ^ J. 



III. ARRAYS OF SU{A) QUANTUM DOTS 

Provided that the conditions outlined above are met, the pillar array of quantum dots may be described by an 
SU(A) invariant Hubbard model. We retain only nearest-neighbor hopping and on-site Coulomb repulsion, and assume 
that no spin-flip or orbital-fiip processes occur. Interdot Coulomb repulsion is not, expected to change our results 
qualitatively. We use one Greek index a — 1, ...,4 to label all four flavors of p-statesEi: \lz = 1,3^ = +1/2) — > |a = 1), 
|1, 1) |2), I — 1, 1) ^ |3), and | — 1, 4) ^ |4). At half-filling the Hubbard Hamiltonian, for an open chain of length 
L sites, may be written: 

H=Y.{^ (4"c,+i,„ + H.c.) + U[n{i) - 2]2 j ; (3) 

here repeated raised and lowered Greek indices are summed over, index i labels the dots, c|" is the creation operator 
for an electron in state a, and n{i) = c\°^Ci „ is the total electron number operator at site i. Both the hopping and 
interaction terms in the Hamiltonian Eq. are explicitly 5t/(4) invariant as can be easily checked by applying a 
unitary transformation, U^cl^ , with WU — 1, which leaves Eq. (^) unchanged. 

At half-filling, the low energy physics of the system is governed by the SU (4) spin degrees of freedom as creation of 
a charge excitation is energetically unfavorable. A weak-coupling renormalization-group (RG) calculation shows that 
umklapp scattering processes (which in the SU (4) case caxry both charge and spin) drive the Hubbard model into a 
Mott-Hubbard insulating phase with gaps in both sectorstj. In the strong-coupling limit of \t\/U <€. 1, again there is 
a charge gap and perturbation theory maps directly the Hubbard Hamiltonian Eq. (^ onto an insulating quantum 
antiferromagnetic Heisenberg spin chain. To 0(t'^/U^) the effective Hamiltonian is: 

^SUW = ^ E I ^°«(^) Tr{Sii)S{i + + 1 MO) [Tr{S{i)Sii + 1)}]4 (4) 

4=1 

plus next-nearest-neighbor terms. Here S^{i) = c\'^Cip — |i5^ are the 15 traceless SU{^) spin generators, the analogs of 
the three Pauli spin matrices in the familiar SU{2) case. With our summation convention Tr{S{i)S{j)} = S'p{i)S^{j). 
The next-nearest-neighbpr terms are 0{t'^ /U^) and tan(6') — C t'^/U'^, where wc find the constant C > 0; its exact 
value csffli be computedta. The purely nearest-neighbor SU{4) spin chain was studied by Affleck, Arovas, Marston and 
RabsoncJ. A combination of exact ground states, RG analysis, and conformal field theory permitted the determination 
of the entire phase diagram, the antiferromagnetic region of which is depicted in Fig. 0. 
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FIG. 2. Antiferromagnetic part of the 5'(7(4) phase diagram for the isotropic, nearest-neighbor 517(4) spin chain. The two 
valence p-electrons on each quantum dot are depicted by filled circles. 5'C/(4) singlet bonds encapsulate four electrons and are 
depicted as rectangles. Chains with an odd number of artificial atoms have free spins at the chain ends in both the dimerized 
and the charge-conjugation (C-breaking) phases. For an even number of atoms, however, the dimerized phase has no free spins. 

There is a spin gap at both weak and strong coupling; we therefore expect the gap to persist at all values of t/U . For 
small values oit/U, 9 is also sipall and in the low temperature limit (T <^ J/kf, 250° K for the In As /Ga As dot) the 
system is in a dimerized phaseEj with broken translational symmetry which can be qualitatively described as a set of 
nearest-neighbor SU (4) singlet bonds as depicted in Fig. ^. Spins not connected by a singlet bond are uncorrelated; in 
other words each of the 6 x 6 = 36 possible configurations of spin and orbital momentum on two such sites are realized 
with equal probability. In contrast, spins on sites connected by a singlet bond are tightly constrained: there is zero 
amplitude for the same configuration to be found simultaneously on both of the sites. This has direct experimental 



consequences as explained below at the end of Sec. IV, The dimerized state, which also breaks reflection symmetry 
about site centers, has a large excitation gap since 0(J) energy is required to break a bond. Consequently spin-spin 
correlations decay exponentially as {Tr{S{i)S{j)}) oc exuL— |i — where ^ is the spin-spin correlation length. 

White's infinite-size density matri:x_renormalization Grouped (DMRG) analysis with open boundary conditions at the 
chain ends confirms this scenarios and determines ^ to be of order the lattice spacing at = 0, see Fig. ^ (a). 
Another quantity of interest here is dimer-dimer correlation function, {Tr{S{i)S{i + l)}Tr{S{j)S{j+l)}), which tells 
us the probability to find a dimer on the link between sites i and i + 1 given that there is one between sites j and 
j -I- 1. The open boundary condition at the chain ends favors one of the two possible dimerization patterns, see Fig. 
H (b). The amplitude of dimer-dimer correlation, the difference between its maximum and minimum values, can be 
used as an order parameter which provides a quantitative measure of the degree of dimerization. 

It is interesting to note that dimer order can be achieved in ordinary translationally-invariant SU (2) antiferromag- 
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netic chains only with large next-nearest-neighbor or biquadratic exchange. Here, however, the dimerized state at 
6* = is a natural consequence of the enlarged SU{A) symmetry. For 6 > 9* = tan^^(l/2) w 0.4636, the chain is in 
a new phase of matter - not realizable for ordinary SU (2) chains - characterized bji spontaneously broken charge- 
conjugation (C-breaking) symmetry, a spin gap, and extended singlet valence bondtcJ. The C-breaking state, unlike 
the dimerized state, breaks reflection symmetry about the centers of bonds, see Fig. ^. We find that the spin-spin 
correlation length increases, and the dimer-dimer order parameter decreases, as the system approaches the transition 
to the C-breaking phase ai = 0*, see Fig. ^ (a) and (b). It may, however, be difficult to reach the C-breaking phase 
in experimental realizations of the system: as t/U is increased, terms in the effective Heisenberg model Eq. (j^ such 
as the nexL-nearest-neighbor exchange^Tr{S'(i)S'(i -I- 2)} become increasingly important. This term has a positive 
coefficientllj, favoring dimerized ordeiEj. 
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FIG. 3. DMRG calculation of (a) the spin-spin and (b) the dimer-dimer correlation functions for odd lattice separation 
and centered at the middle of a spin chain of length L = 36. Diamonds: 8 = 0; circles: 8 = 0.2; and squares: 8 = 0.416. The 
dimerized order diminishes as 8 approaches 6* = tan^^(l/2) ~ 0.4636 



IV. SYMMETRY BREAKING 



It is important to establish whether or not the phases of the pure SU (4) invariant system survive in the presence of 
symmetry breaking processes. We show that the massive dimerized phase is robust in realistic experimental situations. 
The major symmetry breaking process is due to electron-electron interactions which lift the 6-fold degeneracy of the 
configurations of the two p-electrons on each dot. The resulting AU in the Hubbard model induces SU{4) 
SU{2) symmetry breaking both in the on-site energies (6 ^ 3 © 1 © 2) and in the Heisenberg exchange term. Both 
perturbations can be incorporated by perturbing the SU (4) invariant Hamiltonian Eq. (jj) with a two-body interaction 
terms of general bilinear form, which in the simplest case of a translationally-invariant system can be written as: 

L L-l 

H'nund - E Tit; S';{€) + E ^^(*) S'^iz + 1) . (5) 

1=1 1=1 

SU{A) invariance is recovered by setting T^^, T^/f cx (5^(5^; different choices for the tensors then realize all possible 
bilinear SU{4) — > SU(2) symmetry breaking terms. Other less important symmetry breaking processes include: 

1. Non- vanishing hopping between states in neighboring dots with different orbital angular momenta (i^ ^ t S^) 
also breaks S'C/(4) SU{2), as again only spin symmetry remains. The breaking is minimized in the pillar array due 
to rotational symmetry about the vertical axis. 
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2. Spin-orhit, coupling by itself breaks SU{4) SU{2) ® SU{2). However, this effect is small in semiconducting 
quantum dotalj. 

3. Non-magnetic impurities can lift the orbital degeneracy of a dot, breaking SU (4) SU(2), as only spin symmetry 
remains intact. Spin-flip processes, induced by magnetic impurities or external magnetic fields, break SU{-i) all the 
way down to discrete symmetries. It is essential to eliminate both magnetic and non-magnetic impurities in and 
around the semiconducting dots. 

The effects of SU (4) SU (2) symmetry breaking may be analyzed numerically using the DMRG. We find that a 
block size of M = 36 suffices for an accurate description of the massive phases. Even for large values of the symmetry 
breaking parameter corresponding to AU w J the dimer long-range order persists, as is evident in Fig. H. 
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FIG. 4. DMRG calculation of the spin-spin (circles) and the dimer-dimer (squares) correlation functions for odd 
lattice separation and centered at the middle of the chain for the case ^ = and L — 36. We compare the per- 
fectly S'?7(4)-symmetric chain (filled symbols) to one with symmetry broken down to 5(7(2) via AU 7^ (open sym- 
bols). For the broken symmetry case, the on-site tensor T^^ has non-zero entries: T12 = T21 = Tsf — T43 = J/4 and 
rfl = T^^ = Tff = T^i = Tli = Til = Tii = tH = J/4. The nearest-neighbor tensor ff^* has non-zero entries: 

fii = ni = fii = m = J/4 and fit = m = te = m = j/s. 

Other symmetry breaking mechanisms not included in the general bilinear Hamiltonian Eq. (|^) can be incorporated 
by adding a one-body perturbation to the Hamiltonian Eq. (|^): 

L 

1=1 



In particular, spin-orbit coupling corresponds to: 

L 

H'so - A^[5l(^) - Sl{^) - Sl{i) + Sti^)] . (7) 
1=1 

We have examined the effect of this coupling, Eq. (^, and have found that even for an unrealistically large value of 
A = J the dimerization pattern remains intact. The fact that the dimerized phase is robust is not surprising as the 
first excited state is separated by a large, of 0(J), energy gap from the ground state. 

In the extreme limit AU ^ J of large S'C/(4) breaking according to Hund's rules, however, only the triply degenerate 
S — 1 states survive and the chain is describedjia^ an ordinary spin-1 SU{2) quantum antiferromagnet, which is in a 
different massive phase, the Haldanc gap phasec3, with translational symmetry restored as shown in Fig. pi 
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FIG. 5. The Haldane gap phase occurs when the 5(7(4) symmetry is broken down to the usual SU{2) spin symmetry 
of a spin-1 quantum antiferromagnet. SU{2) singlet bonds, depicted as rectangles, involve just two electrons. Translational 
symmetry in the ground state is restored, and there are free spins at the chain ends both for odd and for even chain lengths. 

Transport measurements can be used to confirm the formation of a Mott-Hubbard charge gap at half-filling@0. 
To detect the dimerized spin structure experimentally, it may be possible to exploit the fact that, for an odd number 
of dots only, there are nearly free spins at the chain ends which will dominate the magnetic susceptibility. This 
feature distinguishes the dimerized state from other possible states of the SU{4:) spin chain such as the C-breaking 



ler of sites [see Fig. g]. The free spins may 
ly scanning-tunneling microscopy (STM) 



and Haldane gap phases which have free spins at chain ends for any nuj 
be observable in sensitive electron-spin resonance (ESR) measurementi 
with a magnetized tip, or indirectly via optical spectroscopic experiments 

V. CONCLUSIONS 

Isolated circular semiconducting dots, filled with a few electrons and free of impurities, have already been constructed 
and studiedcl. We propose the construction of a pillar array of such circular dots. Approximate SU (4) symmetry will 
be realized if some simple requirements are met. In particular, the dots must be small to minimize the symmetry 
breaking effect of the intradot electron-electron interaction which partially lifts the degeneracy of the 6 different 
electronic configurations. We predict that a chain of dots, at half-filling (four electrons per dot), will be in an 
insulating, dimerized phase. Four to six dots will suffice because the correlation length is of order the lattice spacing, 
and the state should be robust to various types of symmetry breaking processes as there is a non-zero spin gap to 
low-lying excitations. As a practical application of the proposed quantum dot array there is the problem of quantum 
computationEEl which requires a high degree of quantum coherence between computing elements. The dimerized phase 
is a strongly correlated state and could be used to test the degree of coherence in an array of quantum dots. In 
this it differs greatly from the standard Coulomb blockade seen in coupled dots which operates independently of 
quantum coherence and, apart from the quantization of the electron charge, is classical. Indeed, quantum many-body 
phenomena such as the formation of long-range order are ideal tools to discern quantum coherence. 

Finally we note two possible extensions of this work. Exjperimental evidence for the Kondo effect has been reported 
in transport measurements through a single quantum dotEI. It would be diittiresting to repeat the experiments with 
a SU{4:) dot as increasing the spin degeneracy enhances the Kondo effectEjH. Also, the III-V dots discussed in this 
paper possess enlarged 5'[/(4) symmetry because the 2-fold orbitaLdegeneracy combines with the usual 2-fold spin 
degeneracy. Alternatively, the natural valley degeneracy of silicoro could be esoloited. In Si quantum wells, the 
6- fold valley degeneracy is broken to 2-fold degeneracy by the Si/Si02 interfaceES. This remaining degeneracy, like 
the orbital degeneracy in the III-V dots, is enough to realize overall SU{4) symmetry. 
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